Magnetic relaxation of superconducting quantum dot 
two-dimensional false vacuum decay 

D. R. Gulevich and F. V. Kusmartsev 

Department of Physics, Loughborough University, 
Loughborough, Leicestershire LEll 3TU, United Kingdom 

Abstract 

Quantum tunneling of vortices has been found to be an important novel phenom- 
ena for description of low temperature creep in high temperature superconductors 
(HTSCs). We speculate that quantum tunneling may be also exhibited in mesoscopic 
superconductors due to vortices trapped by the Bean-Livingston barrier. The London 
approximation and method of images is used to estimate the shape of the potential 
well in superconducting HTSC quantum dot. To calculate the escape rate we use the 
instanton technique. We model the vortex by a quantum particle tunneling from a 
two-dimensional ground state under magnetic field applied in the transverse direc- 
tion. The resulting decay rates obtained by the instanton approach and conventional 
WKB are compared revealing complete coincidence with each other. 







1 Introduction 



Quantum tunneling attracts much interest because of its importance for such physical 
systems as helium surface structures, quantum dots, superconductors, that are among 
the most promising for realization of large scale quantum computers. Being one of the 
most important manifestation of quantum mechanics, tunneling process occurs in many 
phenomena in physics, chemistry and biology. 

Large progress has been made in description of tunneling decay when the instanton 
technique was applied. Being a nonperturbative method, it plays a fundamental role in 
description of such type of processes [IJ. It became a powerful tool in physics and has 
got many applications. Some of the advantages of the instanton technique (IT) could 
be: (i) Stronger than standard WKB. Generally, it is hard to tell whether WKB result 
is accurate, whereas the IT is controlled by well defined expansion parameters, (ii) No 
connection formulas, in some cases it is more accurate [2] . (iii) Non-perturbative approach, 
(iv) Instantons, as elementary excitations are topological objects: configurations with dif- 
ferent number of instantons are topologically distinct, (v) Crucial importance for higher 
dimensional field theories. 

The concept of quantum tunneling of vortices in superconductors (e.g. Refs in the 
review ^ about magnetic relaxation in HTSCs) first appeared when measurements of 
magnetic relaxation at ultralow temperatures have been made jlj. The experiments have 
shown that the relaxation rate does not disappear at zero temperature. This phenomenon 
was attributed to the quantum tunneling. Until now it is not well understood yet. Neither 
vortex mass nor Hall coefficient are known exactly. We speculate that the vortex quantum 
tunneling may be also exhibited in mesoscopic superconductors due to trapping potential 
formed by the surface barrier. This would result in higher vortex expulsion magnetic fields 
than expected from pure thermodynamical considerations. 

Due to the Magnus force affecting the vortex dynamics in superconductor, the situation 
is analogous to that of a particle in magnetic field. The problem of a charged particle 
tunneling in presence of magnetic field has itself both theoretical and practical interest. As 
the tunneling is strongly affected by magnetic field, applied in transverse direction, this 
could be efficiently used to control qubits in possible quantum computer realizations of 
the future. In absence of magnetic field the decay rate is related to the imaginary part 
of the free energy as F = \IitiF |_5;. The escape rate can be found making semiclassical 
approximations in the Euclidean path integral. The term corresponding to the ground state 
gives the greatest contribution to the propagator transformed to imaginary times in the 
limit of large time interval. It makes possible to determine the imaginary part of the ground 
state energy. The same considerations must be valid when the magnetic field is applied. 
However, because of the broken time-reversal symmetry, a complex action appears under 
the path integral, when transformed to Euclidean space. Furthermore the imaginary time 
trajectories which extremize the action become complex and the operator corresponding to 
the second variation of the action is a non-hermitian one and as a result possesses complex 
eigenvalues |H] . In this case one could make analytic continuation of the path integral to a 
complex coordinate space or change the time contour in the complex plane. Nevertheless we 
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show that it is possible to make analytic continuation in cyclotron frequency to transform 
the action to a real one. This makes the task much more transparent because of analogy 
with common classical mechanics where everyone has got a good physical intuition. The 
coincidence with the result obtained by usual WKB technique [7| |^ could also serve as a 
proof of validity of the method. Of course it makes sense to make such trick only if we 
consider the task analytically, rather than numerically. Otherwise we could not "return 
back" to real cyclotron frequencies after having got a numerical result with the real action. 

Normally the polar coordinates are used to study the systems that possess rotational 
invariance. However, it becomes hard to work with path integrals in curved coordinates, 
because of additional terms appearing in the action [OlIIIII. Usually one always begins with 
time-sliced path integral in cartesian coordinates before transformed to the curved ones, 
since change of variables in path integrals is not a direct procedure. Everywhere in this 
paper we work with path integrals written in orthogonal coordinates. 

2 Abrikosov vortex in superconducting quantum dot 

2.1 Surface barriers 

The barrier near the surface of type-II superconductors was first studied by Bean and 
Livingstone JT]. It arises from competition of two forces: attraction two the image an- 
tivortex near the border and interaction with the Meissner current. Surface roughness is 
believed to suppress the barrier. However, one can notice that the influence of the surface 
irregularities is much less pronounced for leaving the barrier than for entry 

Another possible source affecting flux dynamics in HTSCs in transverse magnetic field 
is the geometrical barrier |12^ . However, it is not exhibited in superconductors of disk form. 
Moreover, it is expected to dominate magnetic behavior of HTSCs of fiat non-elliptic form 
at elevated temperatures. 

2.2 Dissipation 

There are two main forces affecting vortex dynamics in superconductors: Magnus (Hall) 
force and dissipation. Feigel'man et al. ^H] proposed that the Magnus force is dominant 
in clean superconductors, while other authors argued that the vortex tunneling may 
occur in an intermediate regime. As long as we consider the HTSCs, the dissipative term 
must not be crucial because of small coherence lengths (as well as vortex cores). Indeed, 
the evidence for a low dissipation regime in cuprate superconductors has been presented 
by [ig. 

2.3 Vortex mass 

The same authors [15J argued that the Magnus force is also smaller than standard 
estimates JHl- Thus, the mass of vortex can be relevant to the low-temperature physics of 
clean HTSCs in superclean limit jJT] and should be taken into account in our model. In 
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1965 two contributions to the vortex mass were calculated by Suhl JH]- due to the kinetic 
energy of the vortex core and due to electromagnetic energy. Recently, Chudnovsky and 
Kuklov fTTJ have shown that transversal displacements of the crystal lattice can give a 
significant contribution to the vortex mass. This contribution must be crucial in metals 
with high concentration of superconducting electrons. In our case of small coherence lenght 
^, the most important contribution to the mass arises from the quantization of the electron 
states inside the vortex core (the same paper ^|). It has been shown to exceed the core 
mass by the factor {ep/^Y (Refs [7-9] in IT71). 



2.4 Physical model: summary of basic points 

1. Mesoscopic HTSC disk or quantum dot at low temperature. 

2. Tunneling of a single point vortex trapped by the Bean-Livingston barrier. 

3. No surface roughness: irregularities on the edges are less important for leaving than 
for entry. 

4. No geometrical barrier. 

5. No dissipation: "superclean" limit. 

6. No bulk pinning. 

7. Vortex mass is relevant. 

8. Magnus force is relevant. 

Typical parameters: 

Disk diameter ~ 10 — 100 nm 

Thickness ~ 1 nm 

Coherence length ^ ~ 1 nm 

Bulk penetration depth A ~ 100 nm 



3 The model 

The system of a vortex trapped in HTSC quantum dot is analogous to that of a charged 
particle in 2D potential well with magnetic field applied in transverse direction. For sim- 
plicity we assume h = 1 and m = 1 and consider the limit of large time T (i.e. small 
decay rates), that is usual for the instanton technique. We take the potential in the form 
of rotationally symmetric inverted double well. The case of a particle trapped inside one- 
dimensional inverted double well is studied in details by [12]. See also jH] and PP for ID 
tunneling problems with potential of another shapes. 

The Lagrangian of our 2D model in the Poincare gauge (which coincides with the 
Coulomb gauge in this case) A = (— 0): 

L = — -^{xy -ijx) -U{r) (1) 
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TT/ \ ^2 4 ^ / 2 I 2\ / 2 I 2\2 

U[r) = —r — ar = —[x + y ) — a[x +y) 

Here u denotes the frequency at the bottom of the parabohc potential well, while Uc = 
— = — is the cyclotron frequency. 

The survival amplitude at the bottom of the well expressed in terms of the Feynman 
path integral: 

G(0,T;0,-T) =< 0|e~^^^^|0 >= j Vf{t)e'^-T^^' (2) 
implying the coordinates of the center by 0. Transforming to imaginary times t —>■ —it: 

the Lagrangian transforms to 

L = — V i^{xy - ijx) + f/(r) 

As it was mentioned early the action acquires a complex part after transforming to imagi- 
nary time, in contrast with the case of zero magnetic field. 





Figure 1: Inverted potential - U(r) 



4 Analytic continuation in Uc 

We found useful to make analytic continuation in charge (or uj^. To present some kind 
of grounds to make such a procedure we cite the theorem taken from 120] . 
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Theorem: Consider the next multiple integral: 



F{X,a) 



f{z, a) exp{XS{z, a))dz 



where a = (cti, ...,0;^) is a set of parameters, 7 is a contour in C", with conditions that 
the functions f{z, a) and S{z, a) are analytic, S{z, a) has non-degenerate saddle points at 
Zi, ...,Zs and contour 7 goes through the saddle points when a = Oq. , Then the asymptotics 
of the integral when A — 00 is given by the contribution of the saddle points Zi{a), Zs{a), 



such that Zi[a) 



zi, ...,Zs[a) 



if a is close enough to a^. 



It is straightforward to make generalization of the theorem cited above in the case of 
path integrals, implying "path integral" instead of " n-dimensional integral" and "trajec- 
tory" instead of " n-dimensional stationary point" . We have only one parameter cUc instead 
of the set of parameters a = (ai, ...,0;^) and the condition a = ao corresponds to Uc = 0. 
But we have a real solution of the equation of motion at Uc = 0. It follows from the theorem 
that the asymptotics of the path integral is given by the same analytic formula inside some 
circle around = on the complex plane. Thus if we calculate that one for which the 
stationary trajectories are real ones, the asymptotics inside all the circle can be found by 
continuation in Uc- 

To have even more basis let us look at the convergence of the path integral. Written in 
the sliced form the path integral is: 



Af-l 

dxi...dxN-i I dyi...dyN_i exp — e 

j=0 



ZUJc 



Xj+i - X A / yj+i + yj 




+ U{xj,yj) 



- 2xj+iXj ^ (l/j+i - ^yj+iyj + y]) ^ 



2e 



2e 



lUJr 



+^ [(^J+i - ^j) (%■+! + yj) - (%•+! - i^j+i + ^j)] + {xj, yj) 
Consider the terms with the index k inside the sum: 
^ {2x1 - '^XkXk-\ - 2xk+iXk) + ^ {2yl - 2ykyk-i - 2yk+iyk) + 

iUJc r 

+ [xk+iVk - Xkyk+i - Xkyk - Uk+iXk + UkXk+i + ykXk+ 
+Xkyk + Xkyk-i - Xk-iyk - VkXk - VkXk-i + Vk-iXk] + eU{xk, yk) = 

1 / 2 / \ "^^C / n\ 

= - ( - Xk[Xk-i + Xk+i) - e—Xk{yk+i - Vk-i) j + 

+ ^ fyfe - VkiVk-i + Vk+i) + e-^yk{xk+\ - Xk-iU + eU{xk,yk) 
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It can be seen that the term with ujc must not affect convergence of the Gauss integrals 
because of small factor e in front. 

Thus at least for small we can reduce the stationary point (steepest-descent) method 
to that of Laplace, which is more simple. As a result all the equations are real and the 
instanton trajectories correspond to that ones of a classical particle moving in the inverted 
potential with a transverse magnetic field. 

5 Classical trajectories 

Taking into account the considerations above and making transformation to imaginary 
times t — > —IT as well as uJc — ^ ioJc (in fact the sign in the last procedure must not be the 
matter) we get the following Lagrangian: 




Figure 2: The instanton trajectory slides down from the hill at the center with almost zero 
velocity, bounces from the wall drawing a hint and returns hack to the origin in infinite 
time. There could he many of them, differing hoth in time and angular position. 

The trajectories extremizing the action corresponding to this Lagrangian are that ones, 
which correspond to a classical particle moving in the inverted potential —U{r) and trans- 
verse magnetic field applied in the same direction as before. Now we are going to make 
a little step back from the promise to work in cartesian coordinates and find the classical 
trajectories from the equations of motion written in polar coordinates. Thus: 

L= -^ + —r^(j) + u{r) 
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equations of motion: 

r (0 + —) — const 

We consider the limit T ^ oo as usual for description of ground state decay in the instanton 
approach. The classical trajectories that give the greatest contribution to the path integral 
are that ones, that spend almost all their time at the origin, as the action is zero there. 
This gives E — Q. Let us consider the first equation of motion and the energy conservation 
law: 

r (0 + — ) = const 

+ r^0^ , , 
= ^- \ ^ -U{r) 

Because the trajectory comes from the center of the system, r — > gives const = 0. Suppose 
the opposite. Then at least ~ when r — 0. Obviously this contradicts to the energy 
conservation law as the potential tends to zero under this limit. Hence = — a;c/2 and the 
equations of motion transform to: 

r - n^r + Aar^ = 
with — Lo'^ — ojI/A. The solution of this equation is: 

rci{r) = 



2a cosh Qt 

It corresponds to the instanton with the center at r = 0. There are many other classical 
trajectories with different positions of the centers. Obviously, all of them have the same 
action: 

Wc will take them into account when taking integral over time zero-mode below. Also, there 
are trajectories with the same position in time, but differing from each other by rotation 
around the origin. Similarly, these ones will be counted by the integral over 0-mode below. 
The instanton trajectory transformed to cartesian coordinates reads: 

.,(x) = ..Wc<«(--+^.) = M^ 

Where 1/Q plays a role of the "lifetime" of the instanton, uJc/2 is the frequency of rotation 
around the center. For instance, when lJc/2 ~ fl the instanton makes approximately one 
turn during his "life" , while for bigger uJc/2 the trajectories become spirals spinning around 
the center (although in this considerations cuc is not small already, still it is worthwhile to 
consider such trajectories by the reasons given below). 
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6 Jacobi fields 



Let us calculate the contribution of quantum fluctuations near the classical trajectories. 
Obviously, there is a set of them associated with different 0o- We can fix the one corre- 
sponding to 00 = 0. The other classical trajectories, as well as the fluctuations around 
them, will be taken into account later integrating over the rotation symmetry group. 

In the semiclassical approximation the action is decomposed about the classical trajec- 
tory (for sake of simplicity we omit the normalization constants in front of path integrals): 

J Df(r)e-^[^"1 = e-^=' J D5f{T)e-'^^'^ (6) 

where ^ 

f{T)—rZi{T) + Sf{T) and S'^S — J {Sr,A5r)dT 

The operator A inside the second variation of the Euclidean action is: 



with 



cosh \lT 
+ 41]2 sin2 

^5?nr^ + "^'^"/^ 

U^,, = r5 sm COS = ^ sm(a;cT) 

cosh^Qr 2 2 cosh^^^r ^ ' 

being the second derivatives of the potential evaluated along the classical trajectory. It was 
convenient here to express the frequency cu via and cUc- 

An arbitrary quantum deviation can be decomposed through normalized eigenfunctions 
Xi of the operator A: 

6f{T) = J2CiXi (7) 
Substitution to the path integral leads to the Gaussian integrations over the coefficients Cf 

J J vStt J vStt J vStt 

where det' A denotes product of the eigenvalues of A omitting the zero eigenvalues A^, A,- 
and the negative one A_. They require special treatment and we will pay attention to them 
in the next sections. 

— * — * 

It is convenient to express the resulting survival amplitude (7(0, T; 0, — T) in terms of 
that one for pure parabolic well, that coincides with the contribution of the trivial classical 
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trajectory fli = up to the second order: 

-1/2 



det'A 



det A 



dCr f da 



'It: 



'It: 



1 r-z \ 



(8) 



in this formula Z\ defines single instanton contribution, while is reserved for the trivial 
trajectory. 

In the multidimensional case [22] the ratio of the determinant det' A and det can be 
expressed through that one of the determinants J and Jq of Jacobi fields: 



det'i 



J 



det Aq Jo^(p^T^- 



(9) 



The Jacobi fields satisfy: 



d 



dt \dxidxi 



Jik + 



dxidxi dxidxi 



Jlk + 



d 



dt V dxjdxj 



dxidx] 



Jlk = (10) 



with boundary conditions 



Jik — 0, 



dj 



ik 



Oik 



5. 



ik 



dt m 

The determinant Jq as well as the eigenvalues in the formula ^ will be calculated later. 
Now we begin with evaluation of J. 

The system of four differential equations of the second order (ll(J|) written for the 
Lagrangian (jS)) decouples into two subsystems, each one of the form: 



-i^ - oo^i + u;;r] + u':^^ = 

or in terms of the operator A introduced above: 

A0{r) = 0, with 0{t) = (^^ 

The boundary conditions are: 

e(-T) = 0, i{-T) = l, vi-T) = 0, f,{-T)=0 
for the first subsystem (for which ^ = J^x, i] = Jyx), and 

e(-r) = 0, a-T) = 0, V{-T) = 0, r,i-T) = l 
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(12) 



(13) 



for the second one = Jyx, r] = Jyy). 

Let us find 4 independent solutions of the system. Two solutions of this problem are 
the zero eigenmodes, corresponding to the r- and 0-symmetries. They can be easily found 
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by differentiating the classical trajectory (0) with respect to r and <^o- One can check 
by straightforward substitution that the following solutions are indeed zero-modes of the 
system: 



6-mode: 



r-mode: 



1 /sin ^ 
cosh fir Icos^i^ 



sinh fir 
cosh^ fir 



- cos 
sin 



2 

2 



where we have chosen the time zero mode without the 0-shifting term (obviously it does 
not change anything - any two linear combinations of them could be chosen) and 0o = 0. 

There are two more solutions left. It could seem sophisticated to find them, however 
we will use the following trick. Suppose we deal with a 1-dimensional case. If one solution 
of the second order homogeneous differential equation (written in the Liouville form) is: 



/ 



1 



cosh fir 

then the second one can be found as (for example [SHI or |l()j-2.7.4): 



9 = f 



and if 

the second one is 
9 = f 



dt' 

f = 



dt' 

m 



1 



sinh fir 
2fi 

sinh fir 
cosh^ fir 

sinh^ fir 



+ 



r 



2 cosh fir 



3 sinh fir 
+ -r- 



cosh fir 2fi cosh fir 2 cosh fir fi cosh fir 



The idea to take the anzats in the similar form multiplied by sines or cosines turns out to 
be successful. Indeed, making a straightforward substitution of functions 



sinh fir 



2fi 



+ 



2 cosh fir 



sm 

cos 



2 

UJcT 



sinh fir 3 sinh fir 
2fi cosh fir ^ 2^cosh^ fir 



fi cosh fir 



- cos 
sin 



UJcT 

2 

UJcT 

2 



into the system we conclude that they are indeed the solutions we were looking for. 
Let us analyze the properties of this 4 independent solutions. It is easy to note that: 



'7i(t), V2{'t), ^3{t): ^4(t) are EVEN functions: 
^i(r), C,2{'t), V3{'^)^ V4:{'^) are ODD functions: 



(14) 
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These important properties will be used in future for further calculations. 

At last we are able to find the solution of the system, satisfying the desired boundary 
conditions (fT^ and (fT^ . Expanding 

4 
i=l 

we get the coefficients 



for the vector <^(t) = 




di = -ris{-T), d2 = -r]4{-T), 4 = r]i{-T), = r]2{-T) 

for (f{T) = ( "^^'^ y 

\'^yy/ 

Now let us look at their behavior in the limit T ^ oo. For simplicity we can choose T 
in such a way that sin = and cos = 1 . Obviously, the final result must not depend 
on our specific choice of T, hence choosing T in this way we get significant simplification: 

^i(T) = 0, -2e-^^, 6(T) = 0, ^ --^e^^, (15) 



ryi(T)^2e-^^ V2iT) = 0, r/3(T) ^ ^e^^ V^{T) = (16) 

and 

7^i(-T)^2e-^^, V2i-T) = 0, r/3(-r) ^ -i^e^^ V^i'T) = 
Bearing this in mind we rewrite our coefficients in the following asymptotic form: 

ci = 0, C2~^e^^, C3 = 0, C4~2e-"^ 



Finally, the determinant J can be found: 

J = det 



Jxy Jxy\ 
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7 Elimination of zero eigenvalues 



Now we need to get rid of the zero eigenvalues in the determinant J. This could be done 
by several ways. One could introduce eigenvalue A as a small parameter, perturbating the 
system of differential equations and take the limit A — > of the determinator divided by A^ 
at the end (as we have two zero eigenvalues). Note, that this limit must be taken after the 
limit T oo, as the zero eigenvalues are not exactly zero but tend to it as an exponential 
of T. However this requires the precision at least o(A^). Thus it seems more convenient 
to eliminate A^ and A^ separately using the boundary perturbation method, just looking 
at the behavior of this eigenvalues at large T. As was mentioned above we must expect 
exponential dependence of T. 

Consider the Green function: 

AG{T,T')^-i5{T-T') (17) 

The reason to introduce the minus sign into the definition will be clear below: we will 
receive just the same boundary conditions for the system of differential equations that was 
solved earlier. The boundary conditions for the Green function are: 

G(-r,T')=0, ^(-r,T') = Vr' (18) 
Then the general solution of Aifj = Xqi/j is 

if^lPo- Ao G{t, T')^{T')dT' 

— * A — * 

where ■0o is a solution of the homogeneous equation A'^q — 0. The eigenvalues A^ and At 
can be calculated requiring ip{—T) and ip{T) to be strictly zero (note, that the zero modes 
that we have found early do not satisfy this boundary conditions exactly for finite values 
of T, but only in the limit T — > oo). This leads to the following conditions: 

^/ro(-T)=0 and MT) - J G{T,T')MT')dT' = 

where we have made the Born approximation substituting ip by ipo inside the integral as 
Ao is small. 

First let us find Aq = A^. Then the solution of homogeneous system can be expressed 
via 0- zero mode and other solutions as: 

— * 

and the boundary conditions become 

'^i{-T)+aip2{-T)+(3ips{-T)+'yip,{-T) = 

^,{T) + aip2{T) + (5^,{T) + 7(^1(T) = J^^ G{T, r')^i{r')dr' ^ ' 
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In the last formula the corrections were neglected inside the integral, as they give less 
contribution than (yji (more strictly one can find the values of the coefficients from the 
equations below, substitute them into this integral and prove that it is indeed the case). 
Taking the same T, chosen so that sin = and cos = 1, we get: 

+ P 1 oJ+7 =0 





^ ^ 








' + a \ 


< 




1 \ 








+ a 





)+^Le^T)+7( ]=X^J^^G{T,r'mr')dr' 



Note, that we have the products of two functions of r' inside the integral, that are 
either even or odd ()14|) . Thus only that terms contribute, that consist of functions, both 
even or odd at the same time. Using the explicit form the Green function from the next 
section (|2ip we get the following expression for the integral: 



The two equations from which the desired eigenvalue can be determined are: 

2e-^^ + /3^e^^ = A,^f /^^ {Ur' f + ^,{r'f) dr' 
In the limit of large T this leads to (see Appendix): 

Using the same procedure we find the second zero-eigenvalue A,-. Decomposing the 
solution of the homogeneous system as 

we write the boundary condition: 



(-T) + aifiii-T) + Pifisi-T) + 7V5l(-T) = 



(20) 



Eventually, we get: 

A. ^ 24fi2e-2f^^ 

Note, that A,- has the same form as for the 1-dimensional case (there e — * ^ and the 
operator, for which the eigenvalue is calculated is defined two times smaller than our A), 
but the frequency Vt is dependent on magnetic field now: Vl = u"^ — lo'^ / A. 
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8 Calculation of the Green function 



The Green function, satisfying (fTTjl and p8|l has the next matrix form: 

G(r, r') - 



9u 9i2 

921 922 



where gij = gij{T,T'). 

The equation (fTTjl sphts into two independent subsystems: 



A 

and 



9n\ _ [Sir - t'] 



921 V 



A 



'9i2\ ^ / 
,922) \-6{r - r')^ 

with zero initial conditions following from (fTHj) for both of them. 

Consider the first one. It corresponds to the homogeneous system considered above in 
the regions r < r' and r > r'. In that terms: 9ii{t,t') = (^{r) and 92i{t,t') = rjir). Both 
"left" and "right" solutions must be connected in such a way, that the delta function comes 
out. If we integrate the system over a small interval with the joint inside: 

as ^ and rj are continuous on the joint, the last terms at the left hand sides of the equations 
tend to zero. Obviously, the solution in the "left" region (r < r') is the trivial one: 
^(r)|^<T-' = 0, ri{T)\r<T' = 0. Thus we get the initial conditions for the "right" one: 

'e(r)|.=.'+e = l 

Together with two other conditions coming from the continuity on the border of two regions, 
we have: 

e(r') = 0, ar) = l, Vir')=0, vir') = 

With the replacement r' —T this boundary conditions exactly coincide with ()12p 
(and (jl3p for the second subsystem). Thus the coefficients (that actually are functions 
of r'): 

Cl = -6(^0, C2 = -U{t'), C3 = 6(^')> C4 = 6(7-') 

di = -mi^'), d2 = -rii{T'), ds = ?7i(r'), 4 = 772(^0 
Hence the Green function becomes 

^ ^ \j=l i=l / 
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Using the asymptotic forms ()15|) for t = T we have: 



G(T r') - e(T r') ( '^^''^^^^'''^ + 2e-^^e4(rO -in^^^V^ir') + 2e— v.[r 



40 



4n 
_i_ 

4n 



^"e^^r/i(r')-2e-^'^r/3(r') ^ ^^^^ 



9 Integrals over the peculiar eigenmodes 

Having found the determinator J and the zero eigenvalues, the integrals over zero modes 
left to be evaluated explicitly, as they require special treatment. Recall that the 0- and r- 
modes were found differentiating the classical trajectory fci{T). We can express Tc/ through 
these modes: 



n 



SO that: 



This Teylor series needs to be compared with our expansion ((7j) over the normalized eigen- 
function of the operator A: 



rli + C^X4> + CrXr + ••• = rZi + 



a 



V?2 
\'?2\ 



+ ... 



where the zero eigenmodes were extracted from the sum explicitly. In order to complete 
the integration over the zero modes, the integrals over the coefficients C^, Cr must be 
transformed to the integrals over angular and time positions. The expressions above make 
clear the procedure of changing variables from Ca, Cr to 0, r: 



\Wi\\ 

dCr 
\\V2\\ 



\/2a 



The Jacobian of this transformation is: 



J<h 









ar 


aCr 


aCr 




ar 



det 








a 



where we have have substituted ||<^i|| and ||<^2||5 found to be (see Appendix): 



and 



\V2\ 
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Let us pay attention to the negative eigenvalue. Physical intuition tells that there must 
be only one negative eigenmode, corresponding to the direction in the functional space 
associated with escape of the particle. One can use also oscillation theorems to determine 
how many eigenvalues less than A,- and left. For one-dimensional Shroedinger equation 
this is nothing but number of nodes of the eigenfunction that tells its position number in 
the ordered series of corresponding eigenvalues Aq < Ai < A2 < .... We have managed with 
the negative eigenvalue just as it was done in one-dimensional case (see jT] for details). 
There an analytical continuation of the Gaussian integral over negative eigenmode has 
been used. A special feature is the factor ^, arising from a half of the Gaussian peak. 
The integration over the other half axis turns out to be fake. Thus the contribution of the 
negative eigenmode integral is: 



The sign here depends on how the analytical continuation is done. In future we will omit 
this sign as it is not important. The precise value of A_ is not important also as it cancel 
with A_ in the formula for the ratio of the determinants Q. 
Eventually, 

dCs f dc^ f dc. 1^2, I 1 ,^ r , 

:e~2 = ; ^ / dd) dr 



2ti J V^TT J V^TT 2J-X^2tt y/3a Jo J-t 



We will retain this integrals unsolved for some time. Note, that the integration over the 
angular variable appears naturally from the integration over along the direction in the 
functional space corresponding to the 0-symmetry zero mode. This is nothing to do with 
path integral written in polar coordinates, for which an additional care must be taken. 



10 Contribution of the trivial trajectory 

The determinant Jo of the Jacobi fields along the trivial trajectory fd = is estimated by 
analogy with J. The calculation is rather simple now. 

-UJ.dr -dl + + ^2/4 



An 



So that the system Ao<^(r) = reads: 

-I + uj,r] + + ,^2/4)^ ^ 

with the same boundary conditions (fT^ and (fT!?|l . The system can be solved by standard 
methods. We only write the final answer for the determinants of the Jacobi fields in case 
of the pure parabolic potential: 

Jo 
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11 Instanton gas 



Finally, substituting the determinants and eigenvalues found above to (jHl) and 0, we 
obtain the contribution of one instanton trajectory (more strictly of a group of equivalent 
instantons with different time and angular positions, as we have integrated over the zero 
modes already): 



^0 

with 



Zr 



J 



dr = e'^'^'K / dd) dr 



2^-A_ 27r y/3a Jo J-t 
K 



J~T 



a2TT 

We must also take into account the contribution of the multi-instanton trajectories. Also, 
we must integrate over their positions both in time and angular space. Thus 

_^ = ^-nS,iK^ d<i)i dn d(f)2 dT2... / / rfr„ = e-"^-'ir"^ ^ 



Jo J-T Jo J-T Jo J-T 

represents the contribution of a trajectory with n instantons distributed in the time interval 
{—T,T) and arbitrary directed with respect to the angular coordinate. Finally, we must 
sum all the contributions to get the survival amplitude: 

G(0, T-0,-T) = Zo + Z^ + Z2 + ... = Zof2 e-"^='ir"^^^^^ = Zq exp {27r2TKe~'^^') 

n=0 ^' 

Thus the probability decay rate is: 

F = 2ImE = AirKe-^^' = — e"^^' (22) 

a 

This simple expression is just that one, that we were looking for. To complete the calcu- 
lation we must transform back to the real cyclotron frequencies Uc — ^ —iuJc, so that the 
frequency f2 becomes a square root of the sum now: 

We can rewrite the formula ()22|) in the following way using the expression for the 
classical action (jH): 

F = UnScie-^^' = 12n{^/S^lye-^^' 



The factor v Sd has appeared twice here, just as much as the number of zero modes. This is 
a common situation in the instanton technique (e.g. Coleman, "Aspects of Symmetry" 
p. 337 gets 4 factors ^/Sd in case of four zero modes). 
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12 WKB 

Now we will get the same result using different method. Even the way, in what we 
deal with the magnetic field is different here. The Hamiltonian corresponding to our initial 
system ((H) is: 

2 2 ^ ^' ^ ^ 2 

Using the same gauge A = (— 0), it can be rewritten as follows: 

where L is angular momentum operator: L = xpy — ypx and we have defined 

The important point here is that -f^o introduced above commutes with the angular mo- 
mentum: [Hq, L] = 0. This is because of the rotational invariance of the system, as the 
potential (including the magnetic field term) in the Hamiltonian is a pure function of 
radial distance r = y^x^~+~y^. 

Consider the matrix R{—(f)), rotating the vectors in the counterclockwise direction: 

, cos(f) sin( 

sin (f) cos ( 

Then the angular momentum operator L generates rotations of the wave function: 

This can be seen as follows (e.g. |23 P-219). One can find generator of the rotation unitary 
group by differentiating with respect to 0: 

t^m{-<i^)^u=o = tmR{-<i))f^^R{-m^=o = 

= iVij{f) (^]^^ = -I fyx^ - ^{f) = Lij{f) 

Thus L is indeed the generator of the unitary group of rotations. 
With all this in mind, the amplitude becomes: 

^(0, T; 0, -T) =< 0|e-^2^^|0 >=< Ole''^^"^ e''^^^^ \6 >=< 6\e-'^^°'^\6 > 
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as the rotation of the state |0 > is equal to itself. Thus the problem has reduced to finding 
the decay rate of the particle in the potential of the same shape with the renormalized 
frequency uj — > yjuj"^ + uj"^/ A = VL va. absence of the magnetic field: 



H 



Y^ + Ueff{r) with Ueff{r) = ^-^r'-ar^ 



Shroedinger equation in polar coordinates reads: 



f Id 1 92 



with L = -ih^ and = e^''^/(r): 



' 2 [dr^^ rdr ' + ^-//^^ 



f{r) = Ef{r) 



As we are interested in the tunneling from the lowest state we take / = 0. Consider ip in the 
form /(r) = e^'^^'') with S{r) = Sq — ihSi + ... according to the standard WKB expansion: 

fir) = ei^W 

f{r) = C^S', + S[ + ...)e^nSir) 
f%r)^{-^^S',' + 'is',S[ + y'^ + ..ynSir) 

Substituting to the Shroedinger equation above and equalizing terms near the same power 
of h we get: 

From now we can assume h= 1. The first equation gives: 

So^iJ ^2{Ueff{r)-E)^ J p{r), p{r) = ^2{U,ff{r) - E) 



Dividing the second one by S'q-. 



2S[ + {\nS'Qy + -^0 



that is also easily integrable. Eventually, we get the underbarrier WKB wave function in 
polar coordinates: 

ilj{r) = C = exp I — / p{r)dr 
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where ri is the first classical turning point corresponding to the energy E. The constant 
C can be found comparing the result with the normalized ground state wave function for 
the two-dimensional parabolic well: 

ip[r) = Y — e 2 

For the underbarrier WKB wave function near the bottom (but far enough from the turning 
point, so that WKB works already) the integral under the exponent is approximated by 

p{r)dr ^ J \/02^2 _ 2Edr = — 2Q ~ 2Q ^" E 

for E as a small parameter. Taking E = Q as the lowest energy level for two-dimensional 
parabolic well and equalizing the both wavefunctions, we get: 



^/n^ V TT 



where p{r) was also decomposed up to the lowest order: p(r) ^ Qr. All this leads to 



C 



27re 

The outgoing wave function (r > r2, with r2 as a second turning point outside the well): 
tp{r) = C . exp [ — j p{r)dr + i f p{r)drj 

y/rp{r) V Jri Jr2 J 

This is just the wave function, we are specially interested in. Using it, we can estimate 
density of the probability current: 

id d i C"^ f \ 



' ri 

exp(-2£ PW^r) =^exp(-2£ p(r)dr) 



r 



We must integrate the current density over the circumference to get the total probability 
current coming outside the well: 



D= r%7rrj# = -e-^(^) 
7o e 



As the wave function inside the well is normalized, this expression corresponds to the 
transition probability through the barrier, i.e. decay rate. We have denoted W{E) — 
2 J^^ p{r)dr. To link this result with the one, obtained by the instanton technique, we must 
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represent the decay rate in terms of 1^(0). This can be done as follows. We separate the 
integral into three parts: 



rn rr'2 f-r-i 

W{E)/2= / p{r)dr + / p{r)dr + / p{r)dr 

Jri Jr[ Jr'^ 

Where the points r'-^ and r'2 are chosen between the turning points ri and r2 so, that the 
first integral is calculated keeping the quadratic term only of the potential, the second one 
is calculated making expansion in E/U{r) and the last integral is calculated approximating 
the potential by a linear function near the turning point r2. 

E E . 2r2fi2 



J p{r)dr^ J Vn^r^ - 2Edr = J ^2U,ffir) - — - —\n 



E 



p{r)dr ^ I J2Ueff{r)— / =dr 



eff\ 



2f/e//(r) + |ln 



{n + ^Q? - 2ar'i)r[ 



Near the second turning point r2 we approximate potential by a linear function with zero 
at the the point tq = n/^/2a, such that Ueff{ro) = 0: 

Ueffir) ^ -==(r -ro) 

v2a 

so that the third integral becomes 



2aE . fi3 



e//lOrfr-^J2(-^(ro-r^) 



p{r)dr ^ J^^ y 2 (;^^(^o -r) - E^dr = j^^ 
Combining this terms altogether and taking the limits r'^ — and r'2 Tq we obtain: 

WiE)^WiO)-^-^h/'^^' 



n n \aE 

that gives 

D = l^e-'^^o) 
a 

where 



pro I pro /• 

W{0) = 2 JUeff{r)dr = 2 V^-^r^ - 2ar^dr = 

Jo ^ Jo Jo 



Vfi2 _ 2ar^dr^ = — = 
6a 



what coincides with Sd ©. Evidently, using the different method we have got the same 
result ()22|) that was derived by the instanton technique. 
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13 Discussion 



In deriving the formula for the decay rate using the instanton approach we have re- 
stricted ourselves to the small values of cUc. As the standard WKB results gives the same 
formula, we conclude, that the obtained result is more general. Therefore, it is worthwhile 
to look at the instanton trajectories even for uJc/2 ~ uj. Obviously in this regime the 
trajectories are spirals, spinning out and in near the top of the hill. The question is how 
do the trajectories near the top "feel" the behavior of the potential on large distances and 
reproduce the right dependence of the decay rate on the 4th power term of the double well? 
Let us analyze this situation more carefully. During almost all the period of the instanton 
existence the spirals are perpendicular to the gradient direction whereas the velocities are 
small at the top, because of zero energy at the top. Thus both the trajectories and the 
action must be strongly dependent on high order derivatives of the potential, i.e. the 4th 
power term. 

Although this problem can be reduced to itself under zero magnetic field as was shown 
above, we believe that the calculations aboce could serve at least as an instructive example 
of how the instanton technique works in case of a two-dimensional system with transverse 
magnetic field. 

Thus, (i) The instanton technique can be successfully applied to 2D models in magnetic 
field, (ii) Coincidence with standard WKB was obtained, (iii) Analytical continuation in 
ujc works. 

14 Vortex expulsion problem 

Now we are able to calculate the expulsion magnetic field, at what the vortex leaves 
superconducting dot. Using the London approximation and the method of images, the 
potential well inside the disk can be estimated pi] : 

^ + ln(i?/0 - h{l - r') + ln(l - r^) 

where h = HttR^/^q and r is the distance from the center in units of the radius R. The 
minimum exists until the field h = 1. However, due to quantum tunneling the vortex leaves 
the well at higher h > 1. We approximate this potential by the inverted double well and 
use the formula for the decay rate calculated above to estimate the expulsion field. This 
leads to the next dependance on the radius of the disk on the figure. For comparision 
we also present the typical dependance due to thermal activation process for some fixed 
temperature. 



V{r) 



167r2A2 
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Figure 3: Dependance of the expulsion field on the radius of the HTSC disk 

15 Appendix A: Useful integrals 

dx 2 
yo cosh^ X 

sinh^ X 2 
-dx 



oo cosh^x 3 



°° 1 , 4 

T — dx = - 

oo cosh X 3 

16 Appendix B: Faddeev- Popov procedure 

Zero eigenvalues can be eliminated using the Faddeev-Popov procedure that is very 
useful for systems with constraints in quantum field theory. As our initial path integral ^ 
is not properly defined because of overcounting of zero eigenmodes, we need to define 
correct measure of integration in order not to take into account the contributions from the 
equivalent trajectories related by symmetry group transformations (in the field theory this 
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usually corresponds to gauge transformations). Thus, there are two constraints on our path 
integral: 

Fi[f(r + r/, 00 + e)] =< f(r + r],(f)o + e), ip,{r, 0o) > 
F2[f(r + r]Ao + 0)] =< r(r + r/, 0o + <^2(r, 0o) > 
Consider the identity: 

A[Fi[f(r, 0o)], i^2[r(r, 0o)]] T dd6{F,[r{T + 7^,00 + ^)])<5(i^2[r(r + 7^,00 + ^)]) = 1 



T Jo 



(24) 

where the Faddeev- Popov determinant A[Fi, F2] is the Jacobian of transformation of inte- 
gration variables. A, as well as the action along the trajectories close to the classical ones, 
is invariant under translations in time and 0: 

A[Fi[f(r + r/, 00 + e)lF^[f{T + 7^,00 + 6)]] = A[Fi[r(r, 0o)], F^[r{T, 0o)]] 

Indeed, this can be easily proved by shifting variables inside the integral 
Introducing the identity (j24|l to the path integral (jH)): 



T ^27r 

■T 



dO J Dfe-^["^("''^o)]A[Fi, F2]S{Fi[r{T + r/, 0o + e)])S{F2[f{T + r/, 0o + 9)]) = 
df] r dd [ Z}fe-^[^>''^°)lA[Fi,F2]5(Fi[f(r,0o)])5(F2[f(r,0o)]) 



where we have used the invariance both of the Faddeev-Popov determinant and the action 
as well as the identity Dr(r + r],(j)Q + 6) = Df{T, 0o). Thus we have have transformed our 
integrations over zero modes to the integration over the groups of translational symmetry. 
The Faddeev-Popov determinant is found to be 

/ dFi dFi ^ 

A[Fi,F2] = det^ ^ 



\ ' / 0=0,ri=0 



<T^'^i(^'</'o)> <^,^i(r,0o)>\ _ 



Apf I 900 _ 

<^'^2(r,0o)> <^,^2(r,0o)>; v^, 



a 
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